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Abstract. We review the current status of axions as dark matter. Motivation, 
models, constraints and experimental searches are outlined. The axion remains 
an excellent candidate for the dark matter and future experiments, particularly 
the Axion Dark Matter eXperiment (ADMX), will cover a large fraction of the 
axion parameter space. 
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1. Introduction 

Despite our knowledge of dark matter's properties, what it consists of is still a mystery. 
The standard model of particle physics does not contain a particle that qualifies 
as dark matter. Extensions to the standard model do, however, provide viable 
particle candidates. The axion, the pseudo-Nambu-Goldstone boson of the Peccei- 
Quinn solution to the strong CP problem [U [2l [3l [4] , is a strongly motivated particle 
candidate. We review the strong CP problem and its resulting axion in section O 
In the early universe, cold axion populations arise from vacuum realignment 



UM\I\ and string and wall decay 1 III III [11 [II [H [HI [M [H [H [H . Which 



mechanisms contribute depends on whether the Peccei-Quinn symmetry breaks before 
or after inflation. These cold axions were never in thermal equilibrium with the rest of 
the universe and could provide the missing dark matter. The cosmological production 
of cold axions is outlined in section [3l 

Current constraints on the axion parameter space, from astrophysics, cosmology 
and experiments, are reviewed in section [l] The signal observed in direct detection 
experiments depends on the phase-space distribution of dark matter axions. In 
section [5l we discuss possible structure for dark matter in our galactic halo and touch 
on the implications for detection. Experimental searches and their current status are 
discussed in section [6l Finally, the current status of the QCD axion is summarized in 
section [71 

This is a brief review, designed to overview the current status of axions as dark 
matter. For more advanced details, we refer the reader to the extensive literature (e.g. 
reviews can be found in refs. [2T | [22 l [23 l (24 1 [25] ) . 

2. The strong CP problem and the axion 
2.1. Strong CP problem 

The strong CP problem arises from the non-Abelian nature of the QCD gauge 
symmetry, or colour symmetry. Non-Abelian gauge potentials have disjoint sectors 
that cannot be transformed continuously into one another. Each of these vacuum 
configurations can be labelled by an integer, the topological winding number, n. 
Quantum tunnelling occurs between vacua. Consequently, the gauge invariant QCD 
vacuum state is a superposition of these states, i.e.. 



The angle, 6, is a parameter which describes the QCD vacuum state, 16*). 

In the massless quark limit, QCD possesses a classical chiral symmetry. However, 
this symmetry is not present in the full quantum theory due to the Adler-Bell-Jackiw 
anomaly [26l[2I]- In the full quantum theory, including quark masses, the physics of 
QCD remains unchanged under the following transformations of the quark fields, qi, 
quark masses, m^, and vacuum parameter, 9: 




(1) 



n 




(2) 
(3) 
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(4) 
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where the ai are phases and 75 is the usual product of gamma matrices. This is not 
a symmetry of QCD due to the change in 9. 

The transformations of Eq. ^ through Eq. ^ can be used to move phases 
between the quark masses and 6. However, the quantity, 

6 = 6 — argdet Al = 9 — arg(miTO2. ..mAr) , (5) 

where M is the quark mass matrix, is invariant and thus observable, unlike 6. 

The presence of 9 in QCD violates the discrete symmetries P and CP. However, 
CP violation has not been observed in QCD. An electric dipole moment for the neutron 
is the most easily observed consequence of QCD, or strong, CP violation. 9 results in 
a neutron electric dipole moment of [211 ESI ESI El] , 

|d„| - 10"^^ e cm, (6) 

where e is the electric charge. The current experimental limit is [28] 

|d„| < 6.3 X 10-2^ ecm (7) 

and thus, \9\ < 10^^. There is no natural reason to expect 6 to be this small. CP 
violation occurs in the standard model by allowing the quark masses to be complex 
and thus the natural value of 6 is expected to be of order one. This is the strong 
CP problem, i.e. the question of why the angle 6 should be nearly zero, despite the 
presence of CP violation in the standard model. 

The Peccei-Quinn (PQ) solution [Tl[2] to this problem results in an axion [3 [4]. 
While other solutions to the strong CP problem have been proposed, the presence of 
the axion in the PQ solution makes it the most interesting when searching for the 
dark matter of the universe. Thus, we focus on this solution only in the following. 

2.2. The Peccei-Quinn solution 

The axion is the pseudo-Nambu-Goldstone boson from the Peccei-Quinn solution to 
the strong CP problem [1] [51 [31 H] . In the PQ solution, 9 is promoted from a parameter 
to a dynamical variable. This variable relaxes to the minimum of its potential and 
hence is small. 

To implement the PQ mechanism, a global symmetry, U{1)pq, is introduced. 
This symmetry possesses a colour anomaly and is spontaneously broken. The axion 
is the resulting Nambu-Goldstone boson and its field, a, can be redefined to absorb 
the parameter 6. While initially massless, non-perturbative effects, which make QCD 
6 dependent, also result in a potential for the axion. This potential causes the axion 
to acquire a mass and relax to the CP conserving minimum, solving the strong CP 
problem. 

As there are no degrees of freedom available for the axion in the standard model, 
new fields must be added to realize the PQ solution. In the original, Peccei-Quinn- 
Weinberg-Wilczek (PQWW) axion model, an extra Higgs doublet was used. We review 
this model to demonstrate the PQ mechanism. 

Assume that one of the two Higgs doublets in the model, couples to up- 
type quarks with strengths yf and the other, (j)d, couples to down-type quarks with 
strengths yf, where i gives the variety of up- or down-type quark. We label the up- and 
down- type quarks Ui and di, respectively (rather than g;, as in the previous section). 
With a total of N quarks, there are N/2 up-type quarks and down-type quarks. The 
leptons may acquire mass via Yukawa couplings to either of the Higgs doublets or 
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to a third Higgs doublet. We ignore this complication here and simply examine the 
couplings to quarks. 

The quarks acquire their masses from the expectation values of the neutral 
components of the Higgs, 0^ and 0^. The mass generating couplings are 

- yyl<l>lum + yfdl<l>°ddm + h.c. . (8) 
Peccei and Quinn chose the Higgs potential to be 

V{(f>u,(j)d) = -y.l(j)i4>u-y?d(t>\(t>d^^aij(j)\(j)^(l)](j)j^ ,(9) 

where the matrices (a^) and (6^) are real and symmetric and the sum is over the two 
types of Higgs fields. With this choice of potential, the full Lagrangian has a global 
Upq{1) invariance , 

^ e-*""T=zi, (12) 
->e-*"''T=d, (13) 
e ^6~N{a^ + ad) . (14) 
When the electroweak symmetry breaks, the neutral Higgs components acquire 
vacuum expectation values: 

«)=z;„e'^"/''" (15) 

- t-rfe^^"/"- . (16) 

One linear combination of the Nambu-Goldstone fields, P„ and Pd, is the longitudinal 
component of the Z-boson, 

h = cos [3yPu - sin [3yPd ■ (17) 
The orthogonal combination is the axion field, 

a = sin j3yPu + cos (3yPd ■ (18) 
Using Eqs. p^ . pT]) and pB]) with Eq. ([5]), the axion couplings to quarks arise 

from 

- Cm = m^u\,^e^^^'^URi + mfd^^e* "dffi + h.c. , (19) 

where = y^'-Vu and = yfwd. The axion field dependence can be moved from 
the mass terms using the transformations of Eqs. (|12p . P^ and p^ . Direct couplings 
between the axion and quarks will remain in the Lagrangian, through the quark kinetic 
term. Defining v = ^ v\ + , the corresponding change in Q is 

d^d-N{vulvd^Vdlvu)alv , (20) 
which can be absorbed by a redefinition of the axion field. 

Non-perturbative QCD effects explicitly break the PQ symmetry, but do not 
become important until confinement occurs. These effects give the axion field a 
potential and when significant, the field relaxes to the CP conserving minimum. Hence 
the PQ mechanism, which replaces Q with the dynamical axion field, solves the strong 
CP problem. 

Under the PQWW scheme, the axion mass is tied to the electroweak symmetry 
breaking scale, resulting in a mass of the order of 100 keV. This heavy PQWW axion 
has been ruled out by observation, as discussed in Section l4Tl This does not, however, 
eliminate the possibility of solving the strong CP problem with an axion. In the 
following section, we discuss viable axion models. 
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2.3. Axion models 

"Invisible" axion models, named so for their extremely weak couplings, are still 
possible. In an invisible axion model, the PQ symmetry is decoupled from the 
electroweak scale and is spontaneously broken at a much higher temperature, 
decreasing the axion mass and coupling strength. Two benchmark, invisible axion 
models exist: the Kim-Shifman-Vainshtein-Zakharov (KSVZ) [29l [30] and Dine- 
Fischler-Srednicki-Zhitnitsky (DFSZ) [3ll[32] models. In both these models, an axion 
with permissable mass and couplings results. 

In the KSVZ model, the only Higgs doublet is that of the standard model. The 
axion is introduced as the phase of an additional electroweak singlet scalar field. The 
known quarks cannot directly couple to such a field, as this would lead to unreasonably 
large quark masses. Instead, the scalar is coupled to an additional heavy quark, also 
an electroweak singlet. The axion couplings are then induced by the interactions of 
the heavy quark with the other fields. 

The DFSZ model has two Higgs doublets, as in the PQWW model, and an 
additional electroweak singlet scalar. It is the electroweak singlet which acquires a 
vacuum expectation value at the PQ symmetry breaking scale. The scalar does not 
couple directly to quarks and leptons, but via its interactions with the two Higgs 
doublets. 

PQ symmetries also occur naturally in string theory, via string compactifications, 
and are always broken by some type of instanton. While this could be expected to 
make the axion an outstanding dark matter candidate, string models favour a value of 
the PQ scale that is much higher than that allowed by cosmology (see the discussion 
in [3]) . A review of the current situation can be found in Ref. [33] • As discussed in 
[33] . it is difficult to push the PQ scale far below 1.1 x 10^^ GeV and easier to instead 
increase its value. 

Generically, axion couplings to other particles are inversely proportional to fa, 
however, the exact strength of these couplings are model dependent. For example, the 
coupling between axions and photons can be written, 

^077 — gaijaE ■ B , (21) 

where E and B are the electromagnetic field components. The coupling constant is 

3a77 = y J (22) 

TTJa 

where a is the electromagnetic fine structure constant, 10^ GeV < /a ^ 10^^ GeV is 
the axion decay constant, of the order of the PQ scale, and is a constant containing 
the model dependence. Explicitly, 



2 \N 3(1 

where z is the ratio of the up and down quark masses, N is the axion colour anomaly 
and E, the axion electromagnetic anomaly [34] . The term containing the ratios of light 
quark masses is approximately equal to 1.95. The model dependence arises through 
the ratio, E/N . In grand-unifiable models, E and N are related and E/N = 8/3. The 
DFSZ axion model falls into this category and in this case, ~ 0.36. For a KSVZ 
axion, E = {] and g^ = —0.97. 

It is possible for an axion to solve the strong CP problem, as shown by the 
existence of the KSVZ and DFSZ axion models. While significant for that alone, the 
axion also provides an interesting candidate for the cold dark matter of the universe. 
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3. Cosmological production of axions 

3.1. Properties of axion dark matter 

Axions satisfy the two criteria necessary for cold dark matter: (1) a non-relativistic 
population of axions could be present in our universe in sufficient quantities to provide 
the required dark matter energy density and (2) they are effectively collisionless, i.e., 
the only significant long-range interactions are gravitational. 
Despite having a very small mass [HI [551 [131 [M] , 



axion dark matter is non-relativistic, as cold populations are produced out of 
equilibrium. There are three mechanisms via which cold axions are produced: vacuum 
realignment [3 E [7] , string decay [8l[9l[10l[ni[l2l[l3l[l4l[Ill[l6l[ni[l8]and domain 
wall decay [THl [IHl HO] ■ In this section, we discuss the history of the axion field as the 
universe expands and cools to see how and when axions are produced. Wc also review 
vacuum realignment production in detail, as there will always be a contribution to the 
cold axion populations from this mechanism and, as discussed below, it may provide 
the only contribution. A complete description of the cold axion populations can be 
found in ref. |35j . 

3.2. Topological axion production 

There are two important scales in dark matter axion production. The first is 
the temperature at which the PQ symmetry breaks, TpQ. Which of the three 
mechanisms contribute significantly to the cold axion population depends on whether 
this temperature is greater or less than the inflationary reheating temperature, Tr. 
The second scale is the temperature at which the axion mass, arising from non- 
perturbative QCD effects, becomes significant. At high temperatures, the QCD effects 
are not significant and the axion mass is negligible |36j . The axion mass becomes 
important at a critical time, ti, when rriati ~ 1 |5l[6l[7]. The temperature of the 
universe at ti is Ti ~ 1 GeV. 

The PQ symmetry is unbroken at early times and temperatures greater than Tpg. 
At TpQ , it breaks spontaneously and the axion field, proportional to the phase of the 
complex scalar field acquiring a vacuum expectation value, may have any value. The 
phase varies continously, changing by order one from one horizon to the next. Axion 
strings appear as topological defects. 

If TpQ > Tp, the axion field is homogenized over vast distances and the string 
density is diluted by inflation, to a point where it is extremely unlikely that our visible 
universe contains any axion strings. In the case Tpg < Tp, the axion field is not 
homogenized and strings radiate cold, massless axions until non-perturbative QCD 
effects become significant at temperature, Ti. Agreement has not been reached on 
the expected spectrum of axions from string radiation and there are two possibilities. 
Either strings oscillate many times before they completely decay and axion production 
is strongly peaked around a dominant mode [H [9l [Til [HI O [HI [El [S] or much 
more rapid decay occurs, producing a spectrum inversely proportional to momentum 
[TOl [37] . Rapid decay produces ^ 70 times less axions than slow string decay, leading 
to different cosmological bounds on the axion mass (see section [42]) . 




(24) 
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When the universe cools to Ti, the axion strmgs become the boundaries of N 
domain walls. For iV = 1, the walls rapidly radiate cold axions and decay (domain 
wall decay). If iV > 1, the domain wall problem occurs ^38^ because the vacuum is 
multiply degenerate and there is at least one domain wall per horizon. These walls 
will end up dominating the energy density and cause the universe to expand as 5 oc , 
where S is the scale factor. Although other solutions to the domain wall problem have 
been proposed [18], we assume here that = 1 or Tpq > Tr. 

Thus, if TpQ < Tr, string and wall decay contribute to the axion energy density. 
If Tfl < TpQ, and the axion string density is diluted by inflation, these mechanisms 
do not contribute significantly to the density of cold axions. Then, only vacuum 
realignment will contribute a significant amount. 



3.3. Vacuum realignment mechanism 

Cold axions will be produced by vacuum realignment, independent of Tp. Details 
of this method are discussed below, but the general mechanism is as follows. At 
TpQ, the axion field amplitude may have any value. If Tpq > Tp, homogenization 
will occur due to inflation and the axion field will be single valued over our visible 
universe. Non-perturbative QCD effects cause a potential for the axion field. When 
these effects become significant, the axion field will begin to oscillate in its potential. 
These oscillations do not decay and contribute to the local energy density as non- 
relativistic matter. Thus, a cold axion population results from vacuum realignment, 
regardless of the inflationary reheating temperature. 

To illustrate vacuum realignment, consider a toy axion model with one complex 
scalar fleld, (l){x)^ in addition to the standard model flelds. The potential for 4>{x) in 
our toy model is 

V{<j,)^^m''-vl)\ (25) 
When the universe cools to Tpg ^ Va, (f) acquires a vacuum expectation value, 

{4>) =VaeM^G{x)) . (26) 
The relationship between the axion field, a[x), and 9{x) is 

a{x) = Va9{x) . (27) 
The axion decay constant in this model is 

/a ^ ^ • (28) 

For the following discussion, we set = 1. 

At T A, where A is the confinement scale, non-perturbative QCD effects give 
the axion a mass. An effective potential, 

V{e) ^mliT)f Hi- cos 0) , (29) 

is produced. The axion acquires mass, ma, due to the curvature of the potential 
at the minimum. This mass is temperature, and thus time, dependent due to the 
temperature dependence of the potential [36| . 

In a Friedmann-Robertson- Walker universe, the equation of motion for is 

+ m{t)e - ^V20 + ml{T{t)) sin{0) = , (30) 
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where S{t) is the scale factor and H{t), the Hubble constant, at time t. Near 9 = 0, 
sinie) ~ 6. 

We now restrict the discussion to the zero momentum mode. This is the only mode 
with significant occupation when Tpq > Tr, so the final energy density calculated will 
be for this case. When Tr > Tpg, higher modes will also be occupied. For the zero 
momentum mode, neglecting spatial derivatives, the equation of motion reduces to 

e + SH{t)9 + ml{t)0 = O, (31) 

i.e., the field satisfies the equation for a damped harmonic oscillator with time- 
dependent parameters. At early times, the axion mass is insignificant and 9 is 
approximately constant. When the universe cools to the critical temperature, Ti, 
which we define via 

ma(ri(ii))ti = 1 , (32) 

the field will begin to oscillate in its potential [39] . Given the definition of the critical 
time, ti, in Eg. ((5^ [55]. 

Ti c IGeV (^^^Y^) ■ (34) 

The axion field can realign only as fast as causality permits, thus the momentum of a 
quantum of the axion field is 

Pa{ti) ~ ^ ~ lO-^eV (35) 
ti 

for fa — 10^^ GeV, which corresponds to nia — 6 ^eV, by Eq. (|24|) . Thus, this 
population is non-relativistic or cold. 

The energy density of the scalar field around its potential minimum is 

p=lll9^+mlit)e']. (36) 



and 



^1(0') = ^ ■ (37) 



2 

By the Virial theorem, 

'j2\ „2 iq2\ _ _P_ 

P 
J a 

As axions are non-relativistic and decoupled, 

poc^. (38) 

Thus, the number of axions per comoving volume is conserved, provided the axion 
mass varies adiabatically. 

The initial energy density of the coherent oscillations is 

Pi - lfaml{t,)el (39) 
and 01 is the initial, "misalignment" angle. The energy density in axions today is 

Pa = Pl TTT g3U \ ■ 40 
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Using Eqs. ^ and 



PO = T^fa-;- { ] ^1' (41) 



which imphes the axion energy density. 



^^a- 0.15 I -J—) Ol, (42) 



,10i2GeV, 

using Eqs. (1211), (I33|) and jMl)- 

As the axion coupUngs are very small, these coherent oscillations do not decay 
and make axions a good candidate for the dark matter of the universe. 



4. Constraints on the axion 



Laboratory bounds 

The original PQWW axion would have been of order 100 keV mass, and possessed 
couplings large enough to enable the axion to have been produced and detected in 
conventional laboratory experiments. These included searches for axions emitted 
from reactors, where axions would compete with Ml gamma transitions in radioactive 
decay; nuclear deexcitation experiments; beam-dump experiments; and axion decay 
from 1~ heavy quarkonia states, i.e. J/ip and T. All results were negative, thus 
excluding the original PQWW axion within a decade of its prediction. As these limits 
are much weaker than the current astrophysical upper bounds for both the mass 
and coupling to radiation, the reader is directed to earlier reviews [2TJ [22j |34] , and 
the Particle Data Group Review of Particle Properties for discussion and annotated 
limits gnj. 

Searches for axions have also been performed exploiting the coherent mixing of 
axions with photons in experiments realized with lasers and large superconducting 
dipole magnets. These include both searches for vacuum dichroism and birefringence 
through the production of real or virtual axions in a magnetic field '41], and photon 
regeneration (shining light through walls) [HI |331 El] . Axion- photon mixing will be 
discussed briefly further on, but the all such experiments to date have set limits 
on the axion-photon coupling, gaj-y ~ 2 x 10~^ GeV~^ (weakening significantly for 
rUa > 0.5 meV) orders of magnitude weaker than current astrophysical limits and 
direct searches for solar axions. While there arc no prospects for the polarization 
experiments to compete with these latter {ga-f-y ~ 10~^° GeV~^), a new strategy to 
resonantly-enhance photon regeneration may enable them to improve on these limits 
by up to an order of magnitude 45J , and at least one such experiment is in preparation. 

Torsion-balance techniques have enabled searches for axions through their 
coupling to the nucleon spin, and rigorous bounds have been set on their short-range 
interactions for masses below 1 meV. While impressive tour de force experiments, 
how to relate these limits to expectations from the Peccei-Quinn axion is not 
straightforward [iSj . 



4-2. Cosmological bounds 

Whether by the vacuum realignment mechanism, or by radiation from topological 
strings, the production of very light axions in the early universe implies an increasing 
energy density of the universe in axions for lower masses: fla = Pa/pc oc rria^^^ . As 
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the relative importance of the mechanisms is still controverted, a reliable lower bound 
on the axion mass should obtain where the vacuum realignment contribution becomes 
of order unity, fi^ = 0(1). From Eq. ([^^ . this corresponds to nia ~ 6 /ieV(see section 
13. 3p . although there is uncertainty in this estimate itself, owing to lack of knowledge of 
the initial value of 6 within our horizon; the estimate above based on the assumption 
that this parameter is of order unity. An accidentally small value would drive the 
mass associated with closure density downwards, including arbitrarily small masses. 
On the other hand, within the concordance model, f^DM = 0.23, implying that the 
axion mass should be roughly a factor of 4 higher than cited above. Thus the ADMX 
microwave cavity experiment conservatively began its search campaign at rua ^ 2 /ieV, 
and continues to work upwards. Recent discussions of the cosmological bound from 
vacuum realignment can be found in refs. [471 148j . 

The above vacuum realignment bound applies when inflation has homogenized 
the axion field over our horizon. When the PQ scale breaks after inflation, cold axions 
are additionally produced by string and domain wall decay. These extra contributions 
mean that we require inflation to have produced less axions to avoid overdosing the 
universe and thus, the axion mass lower bound increases (see section l373|) . As the 
spectrum of axions from string radiation is debated, we review the possible bounds. If 
axion strings decay rapidly, giving a spectrum inversely proportional to momentum, 
the lower bound for the axion mass is ^ 15 fxeV [T7]. If the decay is less rapid and 
strings go through many oscillations, an analysis based on local strings [15] gives a 
lower bound of 100 /ieV. A similar analysis based on global strings [16] gives a smaller 
lower bound of 30 /leV. Given that the PQ symmetry is global, it is likely that the 
lower bound of 30 ^eV is applicable if the PQ symmetry breaks after inflation and 
axion strings decay slowly. 

4-3. Astrophysical bounds 

In general, introducing a channel for direct or free-streaming energy loss from a star's 
core accelerates the star's evolution. The core will contract and heat up under the 
influence of gravity when axions (or other exotica) compete with the production 
of strongly trapped photons, whose radiation pressure acts to counterbalance 
gravitational pressure. Furthermore, for each stellar system, axions are excluded only 
over a finite range of couplings. As the axion's coupling is increased in the stellar- 
evolution simulations, the free-streaming lower limit of the axion's excluded couplings 
is reached at a point where deviations from an axion-free model first become noticeable. 
However, as the coupling is further increased, the axions themselves eventually become 
strongly trapped; the upper limit corresponds to the regime where their influence on 
evolution diminishes below the threshold of observation. A comprehensive treatment 
of such constraints on properties of axions and other exotica has been published by 
Raffelt [49j . The two most relevant astrophysical limits in framing the region of interest 
where axions may be the dark matter are described below. 

The most stringent constraint on the axion-photon coupling presently is due to 
horizontal branch (HB) stars, i.e. those that are in their helium-burning phase, within 
globular clusters. Globular clusters provide a cohort population of stars of all the 
same age; those that are seen today are have masses somewhat less than that of 
our sun. By the ratio of the number observed in the HB phase to those in the red 
giant phase, i.e. after exhaustion of core-hydrogen burning but before the helium 
flash, one statistically infers an average HB lifetime. The concordance (within 10%) 
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between the calculated and inferred HB lifetime precludes Primakoff production of 
axions (7 + Ze — s- a + Ze, i.e. axions produced by the interaction of a real plus 
virtual photon) at a level corresponding to an upper bound of ga-yj < 10"^'' GeV~^. 
A definitive study of axion cooling in stars using numerical methods aims to extend 
this analysis [50] . 

The lowest-lying upper bound for the axions mass is due to SN1987a. Axions 
produced by nucleon-nucleon bremsstrahlung {N + N — > A^ + A^ + a) during the 
core-bounce of the protoneutron star would have competed with neutrino emission. 
That the duration of the neutrino pulse observed between the 1MB and Kamioka 
water Cherenkov detectors (19 events over 10 seconds) was in good accord with core- 
collapse models precludes axions in the mass range 10~^ eV < nia < 2 eV. This 
range corresponds to the free-streaming regime; as mentioned previously, for axions 
above this range axions themselves are strongly trapped and thus are not effective in 
quenching the neutrino signal. 

5. Phase-space structure of halo dark matter 

The local velocity distribution of dark matter axions is of vital importance for direct 
detection. The Axion Dark Matter experiment (AD MX), detailed in Section [6?2l and 
ref. [51j . in this volume, uses a microwave cavity to search directly for axions. The 
observed signal is the power output from the cavity due to axion conversion to photons, 
as a function of frequency. The frequency corresponds to the energy distribution of 
axions undergoing conversion. As axions are non-relativistic, the signal frequency is 
given by: 



The local velocity distribution thus determines the signal shape. The signal amplitude 
is determined by the density of axions of a particular energy. Thus, the phase-space 
distribution determines the signal observed. 

We expect that the dark halo of the Milky Way consists of a number of components 
which ADMX is capable of observing: (1) a thermalized component with a Maxwell- 
Boltzmann velocity distribution, (2) discrete flows, from tidal stripping of satellite 
halos or coherent dark matter flows crossing the halo, and (3) overdense regions that 
are not gravitationally bound, known as caustics. 

The ADMX search technique assumes that the rates of change of velocity, velocity 
dispersion and flow density are slow compared to the time scale of the experiment. 
The ADMX medium resolution (MR) channel searches for an isothermal component 
of the halo as dark matter axions. The ADMX high resolution (HR) channel searches 
for signals with a narrow velocity dispersion, such as discrete flows. 

Numerical simulations produce large halos within which hundreds of smaller 
clumps, or subhalos, exist [SH[55]. Tidal disruption of these subhalos leads to flows 
in the form of "tidal tails" or "streams." The Earth may currently be in a stream 
of dark matter from the Sagittarius A dwarf galaxy [54j [55] . This stream of subhalo 
debris satisfles both the requirements of small velocity dispersion and a repeatable 
signal and thus may be detectable by the ADMX HR channel. 

Non-thermalized flows from late infall of dark matter onto the halo have also been 
shown to be expected [56l [57] . The idea behind these flows is that dark matter that 
has only recently fallen into the gravitational potential of the galaxy will have had 




(43) 
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insufficient time to ttiermalize witti the rest of tlie halo and wiU be present in the form 
of discrete flows. There will be one flow of particles falling onto the galaxy for the 
first time, one due to particles falling out of the galaxy's gravitational potential for 
the first time, one from particles falling into the potential for the second time, etc. 
Furthermore, where the gradient of the particle velocity diverges, particles "pile up" 
and form caustics. In the limit of zero fiow velocity dispersion, caustics have infinite 
particle density. The velocity dispersion of cold axions at a time, t, prior to galaxy 
formation is approximately Sva '--^ 3 x lO~^'^{10~^eV/ma){to/t)'^^^ [58], where to is 
the present age of the universe. A flow of dark matter axions will thus have a small 
velocity dispersion, leading to a large, but finite density at a caustic. 

The caustic ring model, under the assumptions of self-similarity and axial 
symmetry, predicts that the Earth is located near a caustic feature [59] . This model, 
fitted to rises in the Milky Way rotation curve and a triangular feature seen in the 
IRAS maps, predicts that the flows falling in and out of the halo for the fifth time 
contain a significant fraction of the local halo density. The predicted densities are 
1.5 X 10~^^ g/cm'^ and 1.5 x 10~^^ g/cm'^ [60], comparable to the local dark matter 
density of 9.2 x 10~^^ g/cm^ predicted by Gates et al. [61]. The flow of the greatest 
density is known as the "Big Flow." 

A general treatment of the phase-space structure of dark matter halos, which does 
not require assumptions of self-similarity or symmetry, has recently been developed 
[62] . This treatment studies the statistics of dark matter caustics in the tidal debris 
remaining from mergers of smaller halos to form galaxies and from the primordial 
coldness of dark matter. While more general than the approach of ref. [60], this 
treatment only results in a statistical distribution and does not give specific predictions 
for our galactic halo. 

Additionally, numerical methods have been developed to study caustics and flows 
in dark matter halos. To date, most numerical simulations are too coarse-grained 
to resolve caustic structure, although its presence can be observed when special 
techniques are used [63] iMl IMl IMl ISl]. The recent work of ref. [68] predicts at 
least 10^ discrete streams near our Sun, although specific predictions are not made 
for the stream densities and velocities. 

It has also recently been shown that dark matter axions can exist in the form of 
a Bose-Einstein condensate (BEG) [69] • If this is the case, the formation of caustics is 
suppressed within the BEG. However, vortices are expected to form at the center of 
galactic halos, due to their net rotation. Within a vortex, axion dark matter will exist 
in the normal phase and fiows and caustics will still be present. 

The possible existence of discrete flows provides an opportunity to increase the 
discovery potential of ADMX. A discrete axion fiow produces a narrow peak in the 
spectrum of microwave photons in the experiment and such a peak can be searched for 
with higher signal-to-noise than the signal from axions in an isothermal halo. If such 
a signal is found, it will provide detailed information on the structure of the Milky 
Way halo. 

6. Experimental searches for dark-matter axions 
6.1. Axion-photon mixing 

As pseudoscalars, axions can be produced by the interaction of two photons, one 
of which can be virtual, 7 + 7* — )■ a , the process being known as the Primakoff 
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effect [70] . This implies that photons and axions may mix in the presence of an 
external electromagnetic field, through the Lagrangian density of Eq. (pij) . In all axion 
searches based on the Primakoff effect to date, E represents the electric field of the 
real photon, and B is an external magnetic field. Although the opposite combination 
is possible, it is vastly easier to produce and support a static magnetic field than the 
equivalent electric field; a 1 T magnetic field being equal to an electric field of 30 
MV/cm in gaussian units. In fact, fields of order 10 T are readily achieved today with 
superconducting magnets. The coherent mixing of axions and photons within magnetic 
fields of large spatial extent enables searches of exceedingly high sensitivity, although 
there is yet no experimental strategy capable of reaching the standard Peccei-Quinn 
axion over the entire open range. 

A general formulism for axion-photon mixing in external magnetic fields, including 
plasma effects, is found in Ref. |7l . 



6.2. The microwave cavity experiment for dark-matter axions 

In 1983, Sikivie proposed two independent schemes to detect the axion based on the 
Primakoff effect [TH [73]. The first was a search for axions constituting halo dark 
matter by their resonant conversion to RF photons in a microwave cavity permeated 
by a strong magnetic field. Tuning the cavity to fulfill the resonant condition, 
hiy = maC^{l + O{0^ ^ 10~^)), and assuming axions saturate the galactic halo, the 
conversion power from an optimized experiment is given by: 

P.,lJ-^, (44) 

where B is the strength of the magnetic field, V the cavity volume and Q is the cavity 
quality factor. The most sensitive microwave cavity experiment (and in fact the only 
one currently in operation) is ADMX at Lawrence Livermore National Laboratory. 
This search has excluded axions of KSVZ axion-photon coupling as the local halo 
dark matter, for a narrow range of masses 1.9 < < 3.4 /leV. 

The anticipated conversion power is miniscule, even for the largest superconduct- 
ing magnets feasible; for ADMX the signal expected is of order 10~^^ watts. Fur- 
thermore, as the experiment will necessitate tuning orders of magnitude of frequency 
in small frequency steps, there are limits on how long one may integrate at each fre- 
quency to improve the signal-to-noise, as governed by the Dicke radiometer equation 



S ^ n ^ Ps^^t ^^^^ 

Ti Pji kTji 

Here s/n is the signal-to-noise ratio, Ai/ the bandwidth of the signal, t the integration 
time, and Pg and P„ the signal and noise power respectively. This puts a clear 
premium on reducing the total system noise temperature, which is the sum of the 
physical temperature and the equivalent electronic noise temperature of the amplifier, 
Tn = Tphys + Teiec- ADMX has recently completed an upgrade from conventional 
heterojunction field-effect transistors (HFETs, or HEMTs) with an noise equivalent 
temperature of Te/ec ~ 2 K, to SQUID amplifiers, whose noise equivalent noise 
temperatures can reach the quantum limit, Teiec ~ 50 mK at 750 MHz when cooled 
to comparable physical temperatures. This strategy will enable them ultimately to 
reach the DFSZ model axions, as well as cover the open mass range much faster. The 
microwave cavity experiments is described by Carosi elsewhere in this volume |51| . 
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6.3. Other current axion searches 

In the same report, Sikivie also outlined how to detect axions free-streaming from 
the Sun's nuclear burning core [751 [73] . Axion production would be dominated by 
the PrimakofF process 7 + ^ a + Ze\ for KSVZ axions, the integrated solar flux 
at the Earth would be given by Fa = 7.4 x 10^^m^[eV] cm~^s~^, emitted with a 
thermal spectrum of mean energy ~ 4.2 keV. For relativistic axions, the conversion 
probability to photons of the same energy in a uniform magnetic field is given by 
P{a — > 7) = n = {\ / A){ga^-^BLY\F{q)\'^ , where B is the strength of the magnetic 
field, and L its lengthfj. F{q) = J dxe^'''^ B{x) / BqL represents the form-factor of the 
magnetic field with respect to the momentum mismatch between the massive axion 
and massless photon of the same energy, q = ka — = (w^ — m^)^/^ — w ~ m^/2cj. 

is unity in the limit ql << 2tt, but oscillates and falls off rapidly for ql > 2tt, 
where the axions are no longer sufficiently relativistic to stay in phase with the photon 
for maximum mixing. 

Utilizing a prototype LHC dipole magnet as the basis for its axion helioscope, 
the CAST collaboration have recently published the best limits on the solar axions, 
Qaj-y < 0.88 X 10"^" GeV~\ valid for rria < 10"^ eV [7S], shghtly more stringent 
than those derived from horizontal branch stars. This collaboration has also pushed 
the sensitivity of the search upward in mass into the region of axion models, by 
introducing a gas (*He) of variable pressure into the magnet bore. In this case, 
the plasma frequency cOp — {iTraNe/me)^^^ = endows the x-ray photon with 
an effective mass; thus full coherence of the axion and photon states can be restored, 
and the theoretical maximum conversion probability achieved for any axion mass, by 
the filling the magnet with a gas of the appropriate density [75] . The mass range can 
thereby be extended upwards in scanning mode, by tuning the gas pressure in small 
steps to as high as feasible. In this manner, axions have now been excluded part- way 
into the Peccei-Quinn model band, gaj-y < 2.2 x 10^^° GeV^^ (95% c.L), valid for 
rua < 0.4 eV [77]. This phase of the experiment continues with '^He gas which will 
permit probing of yet higher masses; for further details, see Zioutas elsewhere in this 
volume [75] . 

Purely laboratory bounds on axions or generalized pseudoscalars have also been 
established without relying on either astrophysical or cosmological sources. In photon 
regeneration ( "shining light through walls" ) axions are coherently produced by shining 
a laser beam through a transverse dipole magnet, and reconverted to real photons in a 
colinear dipole magnet on the other side of an optical barrier [TS] [50] ■ The probability 
to detect a photon per laser photon is given by P{a -^7^0) =11^. While current 
limits from photon regeneration (3077 < 2 x 10"'' GeV~^, nia < 0.5 x 10~^ eV) 
[m 1331 mi have not so far been competitive with solar searches, the scheme may be 
resonantly enhanced utilizing actively locked Fabry-Perot optical cavities to strengthen 
limits potentially by an order of magnitude beyond those of CAST and Horizontal 
Branch stars [45] . 

7. Summary 

The current state of the standard QCD axion is shown in Figure [T] [JD]. While it 
represents an oversimplification of the situation, insofar as the various experimental 

J A useful mnemonic in rate estimates for experiments is that, within a few percent, the factor 
[gwBwLiof ~ 10-1**, where gio = 10'^° GeV-\ Bio = 10 T, and Lio = 10 m. 
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Figure 1. Constraints on the PQ-scale fa and corresponding ma from 
astrophysics, cosmology and laboratory experiments. The light grey regions are 
most model-dependent. From ref. |40] . 



and observational limits on mass are model-dependent, the central point is that there 
is a substantial window for axionic dark matter, and that the upgraded ADMX will 
be able to cover about two of the three decades in mass. One should be open to 
surprise from experiments such as CAST which are looking in regions of mass and 
coupling constant where axions are not expected to be the dark matter, but could find 
axion-like pseudoscalars associated with our first forays into beyond-Standard Model 
phenomenology. 
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